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Chapter  1 
Introduction 

-V 

In  recent  years,  researchers  in  electromagnetics  have  expended  considerable  effort 
capitalizing  on  computational  advances  made  possible  by  new  developments  in  computer 
technology.  These  advances  have  made  it  easier  to  develop  highly  efficient,  specialized 
computer  codes  for  many  scattering  or  radiation  problems.  However,  because  of  the  high 
cost  of  developing  a  code  for  each  specialized  geometry,  it  has  increasingly  become  more 
cost  effective  to  use  possibly  less  efficient,  but  more  general  purpose  codes  which  apply  to 
a  broader  class  of  problems. 

Most  general  purpose  codes  employ  an  integral  equation  formulation  which  is  solved 

c 

by  the  method  of  moments  ^1].  Under  this  type  of  formulation,  the  structure  geometry 
is  usually  modeled  either  as  a  wire  mesh  or  as  a  surface  subdivided  into  discrete  planar 
patches.  The  wire  grid  geometry  modeling  approach  has  been  successfully  used  where 
primarily  far  field  quantities  such  as  RCS  or  radiation  patterns  are  of  interest.  Principal 
advantages  of  the  wire  grid  approach  are  that  the  geometry  is  easily  specified  for  computer 
input,  and  only  one-dimensional  integrals  need  be  numerically  evaluated  in  the  method  of 
moments.  The  wire  grid  modeling  approach,  however,  often  proves  unsatisfactory  where 
near  field  quantities  such  as  surface  currents  or  input  impedance  are  desired.  One  obvious 
difficulty  is  in  interpreting  computed  wire  currents  as  equivalent  surface  currents.  Also, 

,  N, 

however,  the  storage  of  energy  in  the  neighborhood  of  a  wire  mesh  is  not  completely 
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equivalent  to  that  of  a  continuous  surface.  As  a  result,  computed  resonant  frequencies  and 
reactive  components  of  computed  impedances  are  often  shifted  from  their  correct  values. 

Most  of  the  difficulties  involved  in  wire  grid  modeling  may  be  overcome  by  surface  patch 
modeling,  and  severed  such  approaches  have  been  developed.  One  well-tested  scheme  [2] 
makes  use  of  the  electric  field  integral  equation  (EFIE)  formulation  and  a  triangular  surface 
patch  modeling  approach  to  represent  the  geometry.  There  are  two  principal  advantages 
of  this  combination: 

1.  The  EFIE  formulation,  in  contrast  to  the  magnetic  field  integral  equation  (MFIE), 
applies  to  open  structures,  and  it  allows  voltage  and  load  conditions  to  be  easily 
specified  at  terminals  defined  on  the  structure. 

2.  Triangular  patches  are  the  simplest  planar  surfaces  which  can  be  used  to  model 
arbitrary  surfaces  and  boundaries,  and  they  permit  patch  densities  to  be  varied 
locally  so  as  to  model  a  rapidly  varying  geometry  or  current  distribution. 

A  disadvantage  of  the  EFIE  approach,  however,  is  that  the  resulting  integral  equation 
contains  both  a  singular  kernel  and  derivative  operators.  These  features  require  that  an 
algorithm  based  on  the  EFIE  formulation  must  be  developed  carefully  to  ensure  that  it  is 
both  accurate  and  numerically  stable. 

In  this  report,  we  extend  the  approach  of  [2]  to  develop  a  procedure  for  analyzing 
an  arbitrary  configuration  of  conducting  wires  and  bodies.  An  important  feature  of  the 
algorithm  developed  is  its  ability  to  handle  wire-to-wire,  surface-to-surface,  and  wire-to- 
surface  junctions.  Numerous  applications  requiring  these  modeling  capabilities  exist.  For 
example,  many  practical  antennas  are  formed  by  attaching  some  thin  wire  configuration 
to  a  conducting  body,  and  exciting  the  resulting  structure  by  maintaining  a  potential 
difference  between  the  wire  and  conducting  body  at  the  attachment  point.  Often,  the 
shape  or  electrical  size  of  the  conducting  body  is  such  that  it  cannot  be  treated  as  an 


infinite  ground  plane.  Examples  of  such  structures  include  automobile  radio  antennas  and 
wire  antennas  for  shipboard  communications.  In  both  applications,  the  vehicle  supporting 
the  antenna  has  a  complicated  geometry  whose  physical  dimensions  are  only  on  the  order 
of  a  wavelength  at  the  frequency  of  operation. 

Several  approaches  to  specific  wire/body  junction  problems  have  been  reported  in  the 
literature.  Bolle  and  Morganstem  [3]  numerically  solved  the  classically  formulated  case 
of  a  monopole  protruding  from  a  small  sphere  by  considering  it  to  be  the  limiting  case  of 
a  conical  antenna.  Tesche  and  Neureuther[4]  found  current  distributions  on  a  monopole 
mounted  on  a  conducting  sphere  by  using  the  Green’s  function  for  the  electric  field  pro¬ 
duced  by  a  point  source  in  the  presence  of  the  sphere.  Tsai  [5]  solved  the  problem  of  a 
monopole  attached  to  spheres  and  cylinders  based  on  a  Fourier  transform  solution  for  the 
wire  and  the  equivalence  principle  for  the  body.  Cooper  [6]  measured  current  distributions 
and  input  impedances  of  monopoles  attached  to  conducting  cylinders.  Albertsen  et  al.  [7] 
determined  radiation  patterns  for  wires  attached  to  a  circular  cylinder  modeled  by  quadri¬ 
lateral  patches.  Their  formulation  was  based  on  the  EFIE  formulation  for  wires  and  the 
magnetic  field  integral  equation  (MFIE)  for  the  cylinder.  Shaeffer,  Medgyesi-Mitschang 
and  Putnam  [8,9,10]  have  treated  wires  attached  to  bodies  of  revolution  and  to  finite 
cylindrical  bodies.  Richmond  [11]  solved  the  problem  of  a  monopole  antenna  mounted  at 
the  center  of  a  circular  disk  by  using  the  piecewise-sinusoidal  Galerkin  moment  method. 
Marin  and  Catedra  [12]  used  hybrid  moment  method/GTD  techniques  to  treat  a  monopole 
mounted  off-axis  on  a  circular  disk. 

None  of  the  above  procedures  treat  wires  mounted  on  a  surface  at  an  edge  or  vertex. 
Such  problems  seem  to  require  an  a  priori  knowledge  of  the  form  of  the  surface  current 
distribution  near  the  junction.  This  information  usually  must  be  determined  by  obtaining 
a  Green’s  function  associated  with  the  specific  junction  configuration.  Glisson  and  Wilton 
[13,14]  determined  this  distribution  for  an  edge  by  analyzing  a  semi-infinite  current  filament 
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attached  to  a  conducting  wedge.  The  derived  distribution  was  used  together  with  an  EFIE 
formulation  to  treat  a  bent  rectangular  plate  with  an  arbitrarily  oriented  wire  attached 
to  the  plate  at  the  bend.  Newman  and  Pozar  [15]  treated  the  same  problem  by  a  similar 
procedure,  later  extending  it  to  treat  a  wire  attached  to  a  plate  at  or  near  a  knife  edge  [17] 
or  at  a  bend  [18]. 

To  our  knowledge,  none  of  the  above  authors  attempted  to  treat  wires  attached  to  ar¬ 
bitrarily  shaped  vertices  on  conducting  bodies,  presumably  because  Green’s  functions  for 
arbitrary  vertex  geometries  are  unavailable.  Rao  [19],  however,  attempted  to  treat  the  ar¬ 
bitrary  vertex  attachment  problem  by  numerically  determining  the  variation  of  the  current 
near  an  arbitrary  junction  at  an  edge  or  vertex.  However,  his  approach  was  found  to  be 
ill-conditioned  due  to  a  near  linear  dependency  in  the  basis  functions  associated  wiih  the 
wire  junction  and  the  neighboring  surface  patches.  Costa  and  Harrington  [20,21]  closely 
followed  Rao’s  approach,  but  eliminated  the  dependency  in  the  basis  functions.  They  im¬ 
plicitly  assumed  that  for  a  junction  of  any  geometrical  configuration,  only  a  junction  basis 
function  with  a  1/r  variation  is  needed  to  represent  the  surface  current  at  the  attachment 
point.  They  reasoned  that  basis  functions  on  the  surface  would  account  for  any  angular 
variation  of  the  junction  current.  Their  approach  is  appealing  because  of  its  simplicity  and 
generality. 

In  this  work,  we  use  the  approach  of  Costa  and  Harrington — with  a  modification  in  the 
potential  computations  to  increase  the  numerical  efficiency — and  apply  this  procedure  to 
wire-to-surface  junctions  of  arbitrary  geometry.  The  problem  of  modeling  arbitrary  wire- 
to-surface  junctions  is  the  least  well  understood  portion  of  the  present  analysis,  and  hence 
particular  emphasis  is  placed  here  on  validating  numerical  results  from  such  proolems. 
Unfortunately,  bounds  on  computer  resources  limit  our  study  to  structures  whose  largest 
electrical  dimensions  are  only  a  few  wavelengths.  Furthermore,  very  few  measurements  or 
computed  results  involving  antennas  attached  to  edges  or  vertices  may  be  found  in  the 
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literature,  and  hence  we  make  use  of  some  indirect  methods  for  validation  of  the  results. 

In  Chapter  2,  the  formulation  of  the  electric  field  integral  equation  for  an  arbitrary  con¬ 
figuration  of  wires  and  conducting  bodies  is  given.  Current  expansion  functions  and  testing 
procedures  are  introduced  in  Chapter  3  to  convert  the  integral  equation  into  a  matrix  equa¬ 
tion  for  numerical  solution.  Numerical  results  for  a  wide  variety  of  wire- to- surface  junction 
problems  are  presented  in  Chapter  4.  Numerical  examples  are  given  illustrating  the  calcu¬ 
lation  of  current  distribution,  charge  density,  input  impedance,  and  far  field  patterns  for 
various  junction  configurations  including  wires  mounted  on  smooth  surfaces,  edges,  and 
vertices.  A  brief  summary  of  the  research  reported  here  is  found  in  Chapter  5.  Detailed 
derivations  of  the  numerical  methods  used  to  compute  the  potential  integrals  required 
in  the  numerical  procedure  are  given  in  appendices.  A  computer  program,  JUNCTION, 
based  on  the  numerical  procedures  of  this  report  has  been  developed  and  is  described  in 
[22]. 
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Chapter  2 

Electric  Field  Integral  Equation 
Formulation 


In  this  chapter  the  equivalence  principle  is  used  together  with  boundary  conditions  on 
conducting  bodies  and  wires  to  derive  the  electric  field  integral  equation  (EFIE). 

2.1  Equivalence  Principle 

Consider  a  perfectly  conducting  body  which  is  placed  in  a  homogenous  medium  (fi,  e) 
and  whose  surface  is  denoted  by  5  with  unit  normal  vector  n.  The  body  is  illuminated 
by  an  incident  or  impressed  field  (E‘,  H')  due  to  impressed  electric  and  magnetic  current 
sources  J*  and  M\  respectively.  (The  incident  field  ( E\H' )  is  defined  to  be  that  which 
would  exist  in  space  if  the  body  were  not  present.)  The  incident  field  induces  a  surface 
current  on  the  conductor  and  this  current  in  turn  radiates  scattered  fields  (JE*,  H*).  The 
total  field  exterior  to  the  body,  ( E,H ),  is  a  superposition  of  the  incident  field  and  the 
scattered  field  (i.e.,  E  =  E'  +  E*,  H  =  H'  +  H*).  The  induced  surface  current  is  related 
to  the  surface  tangential  magnetic  field  ( J  =  n  x  H),  but  is  unknown.  Our  task  is  to 
derive  an  integral  equation  for  the  induced  surface  current.  If  the  integral  equation  can  be 
solved  for  this  current,  all  other  electromagnetic  quantities  may  be  determined  from  it. 

We  invoke  the  equivalence  principle  [23]  to  allow  the  metal  conductor  to  be  removed  and 
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replaced  by  a  system  of  currents  radiating  in  free  space  such  that  the  fields  produced  are 
identical  to  those  of  the  original  problem.  The  steps  to  generate  this  equivalent  situation 
are: 

1.  First  the  conducting  body  is  removed  and  the  resulting  void  is  filled  with  the  same 
material  as  the  surrounding  medium. 

2.  On  the  mathematical  surface  S  previously  corresponding  to  the  conducting  body, 
we  place  an  electric  surface  current  and  a  magnetic  surface  current  Afeq.  These 
currents  are  required  to  produce  the  same  fields  as  in  the  original  problem;  therefore, 
they  are  called  equivalent  currents,  and  they  are  the  sources  of  the  scattered  field 

3.  The  impressed  currents  as  well  as  the  equivalent  currents  now  radiate  in  an  infinite 
homogeneous  medium  and  the  fields  produced  by  these  currents  can  be  determined 
by  superposition  (or  equivalently,  by  means  of  potential  representations  of  the  fields). 

4.  We  specify  the  exterior  fields  to  be  that  of  the  original  problem  (so  Etma  =  0  on  5). 
We  also  specify  the  field  interior  to  S  to  be  zero. 

5.  The  equivalent  magnetic  current  must  support  any  discontinuity  in  electric  field  at 
S,  but  since  the  tangential  electric  field  is  zero  on  both  sides  of  5,  the  equivalent 
magnetic  current  must  be  zero.  Since  the  interior  magnetic  field  is  zero,  the 
equivalent  electric  current  is  identical  to  the  actual  surface  current  J  on  the  conduc¬ 
tor. 

The  correctness  of  the  equivalent  model  can  be  easily  verified  by  checking  to  see  that  all  the 
postulated  fields  satisfy  Maxwell’s  equations,  including  the  jump  discontinuity  conditions. 
The  above  argument  permits  us  to  assert  that  the  scattered  fields  due  to  currents  on 


conductor  surfaces  can  be  computed  from  potential  integrals  as  if  the  currants  wen  in  an 
infinite  homogeneous  medium.  This  fact  will  be  used  in  the  next  section. 

At  this  point  the  surface  current  J  is  unknown  because  we  do  not  know  the  scattered 
fields  on  the  exterior  of  S.  Since  the  equivalent  current  radiates  in  an  infinite  homogeneous 
medium,  however,  we  may  express  the  scattered  fields  produced  by  J  in  terms  of  potential 
integrals,  and  enforce  boundary  conditions  on  S  to  obtain  an  equation  from  which  J  may 
be  determined. 

2.2  Coupled  Electric  Field  Integral  Equation 

The  previous  section  provides  a  scheme  for  determining  fields  if  the  induced  surface 
currents  are  given.  The  task  now  is  to  derive  an  integral  equation  for  determining  the 
currents. 

Let  S  denote  a  configuration  of  perfectly  conducting  surfaces  immersed  in  an  incident 
electromagnetic  field.  In  general,  S  may  consist  of  a  collection  of  conducting  bodies  and 
wires.  The  wires  may  be  connected  to  smooth  surfaces,  edges,  or  vertices  on  the  bodies, 
forming  different  junction  configurations  as  shown  in  Fig.  2.1.  5  is  then  to  be  viewed  as 
the  union  of  the  bodies  and  the  wires,  which  can  be  expressed  as  follows: 

'  "I 

S  =  Sb  U  Sw  =  (J  SBk  U  jjj  SWi  (2.1) 

where  Nb  and  Nw  are  the  number  of  bodies  and  the  number  of  wires,  respectively.  The 
wire  radii  may  vary  but  are  always  assumed  to  be  small  compared  to  the  wavelength, 
thus  eliminating  the  need  to  consider  any  circumferential  variation  or  components  of  wire 
currents. 

An  electric  field  E\  defined  to  be  the  field  due  to  an  impressed  source  in  the  absence 
of  S,  is  incident  on  and  induces  a  surface  current  J  and  total  current  I  on  Sb  and  Swy 
respectively.  A  pair  of  coupled  integral  equations  for  the  configuration  of  wires  and  bodies 


may  be  derived  by  requiring  the  tangential  component  of  the  electric  field  to  vanish  on 
each  surface.  Thus  we  have 


Etan  =  UUjA  +  V  OH  S 


(2.2) 


where 


A  = 


$  = 


A 

4?r 


r  r  e~>kR  r  //'  e~ikR 

[JsBJ~irdS  +  Jsw2xa(e')  R  dS  ’ 

1  T  t  e~*kR  t  Id/  e.~>kR 

[L  +  JSw 


(2.3) 

(2.4) 


and  R  =  |r  —  r'\  is  the  distance  between  an  arbitrarily  located  observation  point  r  and 
a  source  point  r'  on  5.  In  (2.3)  and  (2.4),  k  =  y,  where  A  is  the  wavelength,  ft  and  e  are 
the  permittivity  and  permeability,  respectively,  of  the  surrounding  medium,  t'  is  the  arc 
length  along  the  wire  axis,  and  a  is  the  radius  of  the  wire.  Eqs.  (2.2)-(2.4)  constitute  the 
EFIE  for  the  unknown  surface  currents  J  and  wire  currents  I.  Because  of  the  thin  wire 
assumption,  we  further  assume  that  only  components  of  electric  field  parallel  to  the  wire 
axis  need  satisfy  (2.2).  In  the  following  chapter,  the  EFIE  is  approximated  by  a  matrix 
equation  for  computer  solution. 
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Chapter  3 

Numerical  Procedures 


In  this  chapter,  a  procedure  for  the  numerical  solution  of  the  electric  field  integral 
equation  is  developed.  Basis  functions  are  chosen  to  represent  the  unknown  currents, 
testing  functions  are  chosen  to  enforce  the  integral  equation,  and  these  are  used  to  derive 
a  matrix  approximant  to  the  integral  equation. 


3.1  Development  of  Basis  Functions 


In  order  to  determine  the  current  distribution  JT,  we  must  first  represent  it  in  a  form 
that  is  convenient  for  numerical  computation.  In  this  section  we  discuss  three  sets  of  basis 
functions  which  may  be  used  to  represent  the  current  induced  on  bodies,  wires,  and  in  the 
neighborhood  of  wire-to-surface  junctions.  It  is  required  that  the  basis  functions  be  linearly 
independent  and  capable  of  approximating  the  actual  surface  current.  A  triangular  patch 
model  of  Sb  and  a  linear  tubular  segment  model  of  Sw  is  assumed  (c.f.  Figs.  3.1  and  3.2). 
Basis  functions  suitable  for  representing  currents  induced  on  Sb  and  Sw  are  given  by 


K(.r)  =  1  hV 


-,  7  =  B  or  W  ,r  in  S 2± 


(3.1) 


0  ,  otherwise 

where,  as  illustrated  in  Fig.  3.3,  S\ J*,  7  =  B  (IF),  is  the  ±  reference  triangle  (segment) 
attached  to  the  nth  non-boundary  edge  (node)  of  a  body  (wire).  The  height  (length) 
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Figure  3.1:  Arbitrary  surface  modeled  by  triangular  patches. 


of  S', J±  relative  to  the  nth  edge  (node)  of  Sy, y  =  B  (W)  is  h J*,  and  p *  is  (±l)x(the 
vector  from  the  free  vertex  (node)  of  S^±  to  r).  The  surface  divergence  of  A^(r),  which  is 
proportional  to  the  linear  charge  density  associated  with  this  basis  function,  is 


V.-^(r)  = 


,  7  =  B  ^mS** 


±ftvv±>  'y =  W  >rinSn± 


(3.2) 


0  ,  otherwise. 

Finally,  we  define  a  basis  function  associated  with  wire  junctions  on  Sb •  A  wire-to- 
surface  junction  is  assumed  to  exist  only  at  a  triangle  vertex.  Referring  to  Fig.  3.4,  the 
vector  basis  function  associated  with  the  nth  junction  is 

(hit)2 


A>)  = 


KrU 

Mr) 


1  - 


(p,+  • 


A&r)  ,rin  SJJ 


m  s;- 


(3.3) 


[  0  ,  otherwise, 

where  the  double  index  nl  refers  to  the  Ith.  triangle,  S^+,  at  the  nth  junction  (c.f.  Fig.  3.4). 
A^i(r)  and  h^+  are  the  previously  defined  body  basis  function  and  the  vector  height, 
respectively,  associated  with  the  edge  opposite  the  junction  vertex  in  S^+.  The  total  flux 
from  the  junction  triangles  into  the  wire  is  normalized  to  unity  if  we  choose 


Knl  = 


Otnl 


Njn 

^ nl  )  ]  ®nl 
l-l 


<Xnl 


(3.4) 


where  a„/  is  the  angle  between  the  two  edges  of  5^+  common  to  the  nth  junction  vertex, 
tni  is  the  length  of  the  edge  opposite  the  junction  in  S^+,  and  aln  is  the  sum  of  the  nth 
junction  vertex  angles.  Nj »  is  the  number  of  patches  attached  to  the  nth  junction.  We 
list  here  some  of  the  properties  which  make  this  basis  function  suited  for  representing  the 
current  in  the  neighborhood  of  a  junction: 
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1.  The  radially  directed  basis  function  has  a  l/|p'|  magnitude  variation  in  the  neigh¬ 
borhood  of  the  junction  vertex.  It  can  be  shown  that  the  current  has  this  radial 
variation  in  the  neighborhood  of  a  thin  wire  junction  regardless  of  the  vertex  shape. 
The  normalization  (3.4)  of  (3.3)  distributes  the  associated  surface  current  approxi¬ 
mately  uniformly  around  the  junction  regardless  of  its  shape;  it  can  be  shown  that 
the  surface  currents  on  Sb  at  the  junction  are  capable  of  representing  any  angular 
variation  of  the  junction  current. 


2. 


The  function  AJn(r )  vanishes  at  the  edge  opposite  to  the  junction  vertex;  thus  no 
charge  is  deposited  along  this  edge. 


3. 


The  surface  divergence  of  AJn(r),  which  is  proportional  to  the  surface  charge  density, 
is  given  by 


V,-^(r)  = 


2Knl 

hi? 


>r  in  sit 


~w  ’rinSi~ 


(3.5) 


4. 


(  0  ,  otherwise. 

The  current  flows  parallel  to  the  two  edges  of  S^f  common  to  the  junction  vertex. 
This  implies  that  no  line  charges  Me  deposited  along  edges  of  S^+. 


The  current  on  the  surfaces  Sb  may  now  be  represented  as 
Nb  Nj 

A')  *  E  +  E  4X(r),  r  on  Sb, 

n=l  n=l 

euid  the  total  axial  current  on  the  wire  may  be  represented  as 


(3.6) 


Nw  Nj 

Hr)t  *  E  C^(r)  +  E  Mr),  r  on  Sw,  (3.7) 

n=l  n=l 

where  Np,  fi  =  B,  W,  or  J  is  the  unknown  number  of  bodies,  wires,  or  junctions,  respec¬ 
tively.  Note  that  according  to  (3.2)  Mid  (3.5),  the  surface  divergence,  and  hence  the  chMge, 
is  constant  on  all  body,  wire,  and  junction  subdomains. 
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3.2  Testing  Procedure 


The  next  step  in  applying  the  method  of  moments  is  to  select  a  suitable  testing 
procedure.  Referring  to  Fig.  3.5,  we  enforce  the  integral  equation  on  Sb  by  integrating  the 
vector  component  of  (2.2)  parallel  to  the  path  from  the  centroid  of  5„+  to  the  middle  of 
the  edge  and  thence  to  the  centroid  of  S®_.  Similarly,  we  enforce  the  integral  equation 
on  Sw  by  integrating  the  vector  component  of  (2.2)  parallel  to  the  path  from  the  centroid 
of  S^+  to  node  m  and  thence  to  the  centroid  of  S%~.  At  a  junction  (c.f.  Fig.  3.6),  we 
first  integrate  the  tangential  electric  field  along  a  path  from  the  centroid  of  each  junction 
triangle  to  the  junction,  and  then  along  the  wire  axis  to  the  center  of  the  attached  wire 
segment.  The  resulting  equations  are  then  combined  into  a  single  equation  for  the  junction 
by  weighting  each  with  the  associated  triangle  vertex  angle,  and  summing  the  results  for 
each  junction  patch.  For  all  path  integrals,  E'  and  A  are  approximated  along  each  portion 
of  the  path  by  their  respective  values  at  the  centroids.  The  integral  on  reduces  to  a 
difference  of  scalar  potentials  at  the  path  endpoints,  thus  the  eliminating  the  requirement 
implied  in  (2.2)  that  $  be  differentiable.  With  l J*,  7  =  B,W,  or  J,  as  the  vector  path 
segment,  we  thus  have 

[  A<+)  ■  C+  +  M<-)  ■  C 1  +  1  *«r)  -  *(0  ] 

=  [  £'«*)  •  C  +  E'(r’-)  •  C  ].  (3.8) 

m  =  1, 2, . . . ,  N^,  7  =  B  or  W, 

and 

-4-  <*ml  )  •  *mt  “  )]  +  j“A(rm  )  *  +  $(**m"  ) 

Qm  /= 1 

“m  l=l 


(3.9) 


m  =  1,2,..  .,Nj. 

Using  the  source  expansion  defined  in  (3.6)  and  (3.7),  we  obtain  from  (3.8)  and  (3.9)  a  set 
of  linear  equations  whose  solution  yields  the  unknown  current  coefficients. 

3.3  Evaluation  of  Matrix  Elements 

Substituting  the  current  representations  (3.6)  and  (3.7)  into  (3.8)  and  (3.9),  we  obtain 
an  N  x  N  system  of  linear  equations,  where  N  =  Nb  +  Nw  +  Nj,  which  may  be  written 
in  matrix  form  as 


■  [zBB]  [zBW]  [zBJ]  ■ 

■  M  ■ 

■  M ' 

[zWB]  [zww]  [zWJ] 

M 

= 

M 

[zJB\  [zm]  [zJJ] 

■  M . 

.  M . 

where  the  elements  of  the  submatrices  are  given  by 

zzi  =  >«[42f-C+42r-«-)+  (3.u) 

1  Njtn 

zii  =  -r  E  <*«i  -  Ku)  +  (3.12) 

/=1 


=  «*:?), 

*!?  = 

♦XT  =  *!(•*), 

(3.13) 

ds' 

(3.14) 

*»r>  -  -)4 Lfs 

t-jkR 

f  V.  ■  AHr'f-g-iS' 

(3.15) 
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7  =  B,W,  or  J  0  =  B,W,  or  J  (3.16) 

and 

vz  —  [£‘(o-c  +  tfwn-cl  (3-17) 

1  Njm 

Vi  =  -J-  £  Qm, £'>£).  +  E:(ri~)-rm-.  (3.18) 

am  1=1 

For  plane  wave  scattering  problems,  for  example,  we  may  set  (c.f.  Fig.  3.7) 

£?(r)  =  (Efi  +  £;d>V'k-r  (3.19) 

where  the  propagation  vector  k  is 

k  =  &(sin  6 ’  cos  <f>'x  +  sin  6'  sin  <j>'y  +  cos &z)  (3.20) 


A  |  A  | 

and  (6  ,<f>)  define  the  angle  of  arrival  in  the  usual  spherical  coordinate  convention,  and 
can  be  expressed  as 

d'  =  cos  8'  cos  <f>'x  +  cos  8i  sin  <f>'y  —  sin  6'z 
<j>  =  —  sin  <f>'x  +  cos  <f>'y. 

For  antenna  radiation  problems  E'e  =  E+  =  0,  in  this  case  the  wire  segments  attached 
to  node  (junction)  m  may  be  thought  of  as  separated  by  a  gap  across  which  the  voltage 
V™  (V^)  is  specified. 

Solution  of  the  linear  system  of  equations  (3.10)  yields  the  set  of  unknown  current 
coefficients  used  in  the  representation  of  the  surface,  wire,  and  junction  currents,  (3.6)  and 
(3.7).  Once  the  surface  current  is  known,  the  scattered  field  or  any  other  electromagnetic 
quantity  of  interest  may  be  readily  determined. 
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Chapter  4 

Numerical  Results 


In  this  chapter  numerical  results  such  as  the  surface  current  near  a  junction  or  the 
input  impedance  of  a  monopole  or  far  field  pattern  are  presented  for  different  junction 
configurations.  The  geometries  considered  include  wires  attached  to  a  smooth  surface,  an 
edge,  and  a  vertex. 

4.1  Calculation  of  Current  and  Charge  Distributions 

Fig.  4.1  illustrates  the  current  distribution  near  the  junction  between  a  circular  cone 
and  a  monopole  attached  at  the  vertex  and  inclined  at  an  angle  of  60°  to  the  axis.  The 
length  of  the  monopole,  the  height  of  the  cone,  and  the  diameter  of  the  cone  base  are 
all  a  =  0.333A.  The  radius  of  the  monopole  is  r  =  0.001a.  The  result  is  compared  with 
the  magnetostatic  current  distribution  on  an  infinite  cone  with  a  semi-infinite  filamentary 
current  attached  to  its  vertex  [24].  The  current  filament  is  also  inclined  at  an  angle  of  60°  to 
the  cone  axis.  Since  the  magnetostatic  result  gives  the  shape  but  not  the  magnitude  of  the 
current  distribution,  a  complex  normalizing  constant  is  chosen  so  that  the  magnetostatic 
current  interpolates  the  dynamic  result  at  one  point.  Also  shown  is  the  uniform  distribution 
that  would  result  from  taking  the  current  in  the  wire  and  distributing  it  uniformly  about 
the  cone. 
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Fig.  4.2  illustrates  similar  quantities  near  the  junction  of  a  circular  disk  and  a  monopole 
attached  at  the  center  and  inclined  at  an  angle  of  30°  to  the  disk  axis.  The  length  of  the 
monopole,  and  the  radius  of  the  disk  are  a  —  0.333A.  The  radius  of  the  monopole  is 
r  =  0.001a. 

Figs.  4.3  and  4.4  show  the  computed  current  distributions  on  a  circular  cylinder  with 
monopoles  attached  to  the  center  of  each  endcap  and  driven  at  the  attachment  point. 
Measured  results  by  Cooper  [6]  are  shown  for  comparison.  The  diameter  and  the  height 
of  the  cylinder  are  0.25A.  The  wire  radius  is  0.007A  and  the  lengths  are  0.125A  and 
0.25A  in  the  two  figures,  respectively.  Similar  results  are  shown  in  Fig.  4.5  for  a  cylinder 
whose  diameter  and  height  are  0.5 A  with  a  wire  length  of  0.125A.  Charge  distributions 
corresponding  to  these  three  cases  are  shown  in  Figs.  4.6,  4.7  and  4.8. 
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4.2  Calculation  of  Input  Impedance 

4.2.1  Wire  mounted  on  a  smooth  surface 

Fig.  4.9  shows  the  input  admittance  as  a  function  of  frequency  for  an  monopole  at¬ 
tached  to  the  center  of  a  flat  plate  and  fed  at  the  attachment  point.  The  length  of  the 
monopole  is  0.421  m,  its  radius  is  0.0008  m,  and  it  is  oriented  normal  to  the  0.914  m 
square  plate.  The  results  are  compared  with  those  calculated  and  measured  by  Newman 
and  Pozar  [15]. 

In  Fig.  4.10  is  shown  the  input  impedance  of  a  monopole  of  length  h  =  d  and  radius 
a  =  0.0165d  mounted  on  a  sphere  of  radius  d  as  a  function  of  the  frequency.  The  results 
are  compared  with  calculations  by  Bolle  and  Morganstem  [3],  [4]. 

In  Fig.  4.11  is  illustrated  the  input  impedance  versus  frequency  of  a  monopole  inclined 
at  an  angle  of  60°  from  normal  and  attached  off-axis  on  a  circular  disk.  The  results  are 
compared  with  measurements  by  Marin  and  Catedra[12].  The  agreement  is  very  good 
except  for  the  reactance  at  the  highest  frequency  measured. 

Fig.  4.12  illustrates  the  input  admittance  of  a  monopole  mounted  near  the  corner  of  a 
box  on  an  infinite  ground  plane.  The  box  is  a  100  mm  cube  and  the  monopole  length  is 
60  mm  with  a  radius  of  0.8  mm.  The  measured  results  of  [25]  on  a  ground  plane  of  780  mm 
x  780  mm  are  also  shown  in  the  figure.  Also  shown  are  results  computed  by  replacing  the 
monopole  with  an  equivalent  tape  model  and  then  modeling  the  structure  as  in  [2]. 


Figure  4.1:  Current  distribution  near  the  junction  of  a  circular  cone  and  a  monopole 
attached  at  the  vertex  and  inclined  at  an  angle  of  60°  to  the  axis.  The  length  of  the 
monopole,  the  height  of  the  cone  and  the  diameter  of  the  cone  base  .are  all  n  =  0.333A. 


Current  distribution  near  the  junction  of  a  circular  disk  and  a  monopole 
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Figure  4.3:  Current  distribution  for  a  dipole  whose  arms  are  attached  at  the  center  of 
the  endcaps  of  a  circular  cylinder.  The  wire  radius  and  length  are  0.007A  and  0.125A, 
respectively. 


Figure  4.5:  Current  distribution  for  a  dipole  whose  arms  are  attached  at  the  center  of 
the  endcaps  of  a  circular  cylinder.  The  wire  radius  and  length  are  0.007A  and  0.125A, 
respectively. 


distribution  for  a  dipole  whose  arms  are  attached  at  the  center  of 
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Figure  4.7:  Charge  distribution  for  a  dipole  whose  arms  arc  attached  at  the  center  of 
the  cndcaps  of  a  circular  cylinder.  The  wire  radius  and  length  are  0.007A  and  0.25 A, 
respectively. 


Figure  4.8:  Charge  distribution  for  a  dipole  whose  arms  are  attached  at  the  center  of 
the  endcaps  of  a  circular  cylinder.  The  wire  radius  and  length  are  0.007A  and  0.125A, 
respectively. 


Figure  4.9:  Input  admittance  as  a  function  of  frequency  for  an  monopole  attached  to  the 
center  of  a  fiat  plate  and  fed  at  the  attachment  point.  The  length  of  the  monopole  is 
0.421m,  its  radius  is  0.0008  m,  and  it  is  oriented  normal  to  the  0.914  m  square  plate. 


Figure  4.10:  Input  impedance  of  a  monopole  of  length  h  =  d  and  radius  a  =  0.0165d 
mounted  on  a  sphere  of  radius  d  as  a  function  of  the  frequency. 


attached  off-axis  on  a  circular  disk. 


Figure  4.12:  Input  admittance  of  a  monopole  mounted  near  the  corner 
plane.  The  box  is  a  100  mm  cube  and  the  monopole  length  and  radiu 
mm,  respectively. 


4.2.2  Wire  mounted  on  an  edge 

Fig.  4.13  shows  the  input  admittance  as  a  function  of  bend  angle  for  a  monopole  attached 
to  the  edge  of  a  plate  and  fed  at  the  attachment  point.  The  length  of  the  monopole 
is  0.25A,  its  radius  is  0.001A,  and  it  is  mounted  at  the  center  of  the  longer  edge  of  a 
0.4A  x  0.5 A  rectangular  plate.  The  reactance  is  not  in  agreement  with  [18],  in  which 
both  an  exact  Green's  function  to  model  the  current  in  the  neighborhood  of  the  junction 
and  a  rectangular  patch  model  was  used,  but  several  different  calculations  have  been  used 
to  validate  the  present  results.  First,  as  shown  in  Fig.  4.14,  the  magnetostatic  current 
distribution  at  the  junction  was  calculated  and  used  to  deduce  weighting  coefficients  for 
the  junction  triangle  basis  functions  for  representing  the  angular  variation  of  the  junction 
current.  However,  this  produced  an  insignificant  change  in  the  result.  Secondly,  the  wire 
was  replaced  by  equivalent  tape  models,  shown  in  Fig.  4.15,  in  which  the  width  of  the 
tape  is  four  times  the  radius  of  the  wire  antenna  [26,  p.  20],  and  either  1,  2  or  3  patches 
are  used  at  the  junction  to  connect  the  tape  to  the  surface.  As  Fig.  4.14  shows,  as  the 
number  of  patches  increases,  allowing  the  current  to  flow  more  freely  in  any  direction  on 
the  plate,  the  results  approach  those  of  the  present  method.  Finally,  we  also  observe  that 
as  the  wire  angle  a  approaches  zero,  the  wire  and  ground  plane  configuration  approaches  a 
quarter  wavelength  section  of  open  circuited  transmission  line,  whose  input  resistance  and 
reactance  should  approach  zero.  The  present  results,  contrasting  to  those  of  [18],  more 
nearly  follow  this  trend. 

4.2.3  Wire  mounted  near  or  at  a  vertex 

Fig.  4.16  illustrates  the  variation  of  the  impedance  as  a  function  of  attachment  position  for 
a  quarter  wave  monopole  mounted  near  or  at  a  corner  of  a  0.4 A  square  plate.  In  Fig.  4.17, 
the  same  monopole  is  attached  at  or  near  a  corner  formed  by  three  0.4A  square  plates.  The 
antenna  is  attached  normal  to  the  top  plate  and  driven  at  the  attachment  point.  Its  radius 
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is  a  =0.0015A.  The  results  for  these  two  cases  are  compared  with  calculated  and  measured 
results  by  Newman  and  Pozar  [15].  The  agreement  is  good  in  all  three  cases.  Since  the 
present  method,  in  contrast  to  [15],  should  also  apply  for  the  wire  located  directly  at  the 
plate  corners,  we  give  results  for  these  cases  also. 

4.3  Calculation  of  Far  Field  Patterns 

Fig.  4.18  shows  the  radiation  pattern  of  a  wire  of  length  0.25A  and  radius  0.007A 
attached  to  a  circular  cylinder  with  diameter  and  height  0.25A.  Figs.  4.19-4.21  show 
similar  results  for  cylinders  of  various  sizes.  The  measured  results  shown  for  comparison 
are  by  Cooper  [6].  In  Fig.  4.22  is  depicted  a  0.25A  monopole  mounted  near  the  corner  of  a 
square  plate  with  1.0 A  edge  length.  The  results  of  Pozar  [18]  sue  also  shown  for  comparison. 
Figs.  4.23  and  4.24  show  a  0.25A  monopole  mounted  on  a  sphere  with  radius  0.2A  and  0.3A 
respectively.  The  results  are  in  good  agreement  with  those  of  Tesche  [4]. 
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Figure  4.17:  Input  impedance  as  a  function  of  attachment  position  for  a  quarter  wave 
monopole  mounted  near  or  at  a  corner  formed  by  three  0.4A  square  plates.  The  antenna 
is  attached  normal  to  the  top  plate  and  driven  at  the  attachment  point.  Its  radius  is 
a=0.0015A. 


Figure  4.18:  Radiation  field  pattern  for  a  dipole  whose  arms  are  attached  at  the  center 
of  the  endcaps  of  a  circular  cylinder.  The  wire  radius  and  length  are  0.007A  and  0.25A, 
respectively. 
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"lgure  4.19:  Radiation  field  pattern  for  a  dipole  whose  arms  are  attached  at  the  center 
>f  the  endcaps  of  a  circular  cylinder.  The  wire  radius  and  length  are  0.007A  and  0.5A, 
espectively. 


Figure  4.20:  Radiation  field  pattern  for  a  dipole  whose  arms  are  attached  at  the  center 
of  the  endcaps  of  a  circular  cylinder.  The  wire  radius  and  length  are  0.007A  and  0.25A, 
respectively. 


Figure  4.21:  Radiation  field  pattern  for  a  dipole  whose  arms  are  attached  at  the  center 
of  the  endcaps  of  a  circular  cylinder.  The  wire  radius  and  length  are  0.007A  and  0.25A, 
respectively. 


Radiation  field  pattern  for  a  0.25A  monopole  mounted  on  a  sphere  of  radius 


Radiation  field  pattern  for  a  0.25A  monopole  mounted  on  a  sphere  of  radius 


Chapter  5 
Summary 


The  electric  field  integral  equation  (EFIE)  is  used  with  the  moment  method  to  develop 
a  simple  and  efficient  numerical  procedure  for  treating  problems  of  scattering  and  radi¬ 
ation  involving  arbitrarily  configured  and  arbitrarily  connected  bodies  and  wires.  Three 
triangular-type  basis  functions  are  used  to  represent  the  physical  currents  on  bodies,  wires, 
and  wire-to-body  junctions.  A  junction  basis  function  is  developed  which  appears  to  be 
applicable  to  any  junction  configuration.  Numerical  results,  which  include  comparisons 
with  calculations  and  measurements  found  in  the  literature,  illustrate  the  versatility,  ac¬ 
curacy,  and  efficiency  of  the  method  for  various  junction  configurations,  including  wires 
mounted  on  smooth  surfaces,  at  edges,  or  at  vertices. 


Appendix  A 

Computation  of  Potential  Integrals 
for  Wire  Basis  Functions 


The  basis  function  associated  with  current  crossing  the  nth  node  of  a  segmented  wire 
model  is  given  by  (3.1),  and  its  divergence  is  given  by  (3.2).  Since  dS'  =  a  d<j>  dt\  the 
vector  and  scalar  partial  potentials  due  to  a  unit  current  source  associated  with  the  nth 
node  are  given  in  terms  of  (3.14)  and  (3.15)  as  (c.f.  Fig.  3.3(a)) 
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where  t  is  arc  length  measured  from  a  segment  endpoint  along  the  axis  of  the  wire  segment 
in  the  positive  current  reference  direction,  and 
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is  the  so-called  exact  kernel  in  which  R  is  the  distance  between  observation  point  r  and 
source  point  r'  on  the  wire  surface. 

A.l  Integrals  Associated  with  a  Segment 

It  is  computationally  efficient  to  simultaneously  compute  all  potential  integrals  associated 
with  a  given  segment  Sw*  (either  S™+  or  S™~).  The  vector  and  scalar  potentials  associated 
with  a  segment  involve  the  following  three  integrals: 

I,  =  •'  K(s,s')  is',  (A.4) 

h  =  (h  ~  s')K(s,s')  is',  (A.5) 

I3  =  t  K(s,  s')  is',  (A.6) 

Jo 

where  s'  is  a  local  segment  arc  length  measured  from  the  end  of  the  axis  of  the  source 
segment,  and  h  is  the  length  of  the  segment.  The  first  and  second  integrals  above  are 
involved  in  the  vector  potential  computation  and  the  third  is  related  to  the  scalar  potential 
calculation. 

A. 2  Treatment  of  Singularities  in  Wire  Potential  In¬ 
tegrals 

In  numerical  evaluation  of  (A.4)-(A.6),  two  cases  are  considered: 

(a)  The  observation  point  is  at  least  two  segment  lengths  away  from  the  source 
segment. 

(b)  The  observation  point  is  less  than  two  segment  lengths  away  from,  or  possibly 
on,  the  source  segment. 

For  case  (a),  the  so-called  reduced  kernel, 
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Kr(s,s')  = 


e-jkRr 

Rr  ’ 


(A.7) 


where 


Rr  =  y/\  r-rf +  a> 


(A.8) 


and  a  is  the  source  segment  radius,  may  be  used  in  place  of  (A.3)  to  eliminate  the  numerical 
integration  which  computation  of  the  latter  implies.  With  this  replacement,  (A.4)-(A.6) 
may  be  numerically  integrated  by  Gaussian  quadrature. 

For  case  (b),  the  kernel  is  rewritten  as 


1  F  r*  e~’kR  —  1  r*  1 

K^  =  r4L—R-d*  +  J.Ad 4 


(A.9) 


The  first  integral  of  (A.9)  now  has  a  slowly  varying  integrand  in  which  the  reduced  kernel 
approximation  can  be  used  (i.e.,  R  is  replaced  by  Rr  of  (A.8))  to  obviate  the  need  for 
double  integration: 


1  r*  e-ikR  _  x 
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The  integral  resulting  from  substituting  (A.10)  for  this  slowly  varying  part  of  the  kernel 
in  (A.4)-(A.6)  can  be  easily  evaluated  numerically  by  Gaussian  quadrature. 

The  second  integral  in  (A.9)  contains  the  rapidly  varying  part  of  the  kernel  for  ob¬ 
servation  points  close  to  or  on  the  source  segment.  Often  the  reduced  kernel  is  used  for 
this  term  as  well,  and  its  contribution  can  then  be  determined  analytically.  For  segment 
lengths  smaller  than  about  five  times  the  wire  radius,  however,  this  approach  is  not  suffi¬ 
ciently  accurate.  On  the  other  hand,  it  is  easy  to  show  that  the  second  integral  in  (A.9) 
is  expressible  in  terms  of  the  complete  elliptic  integral  of  the  first  kind.  To  demonstrate 
this,  we  define  a  number  of  geometric  parameters  associated  with  the  segment  and  the 
observation  point  in  Fig.  A.l.  The  distance  between  observation  point  r  and  source  point 
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R  =  |r  -  r'|  =  J  (sm  -s')2  +  p2  +  a2  -  2pacos<p, 


(A.11) 


where 
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Vector  *  is  the  unit  vector  pointing  in  the  direction  of  the  source  segment,  and  am  is 
the  radius  of  the  observation  segment.  Derivation  of  the  relationship  between  the  singular 
integral  of  (A. 9)  and  the  complete  elliptic  integral  of  the  first  kind  may  now  be  summarized 


as  follows: 
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and  K(/3)  is  the  complete  elliptic  integral  of  the  first  kind,  which  can  be  computed  from 
(17.3.34)  of  [27]. 


The  elliptic  integral  becomes  logarithmically  singular  when  the  observation  point  ap¬ 
proaches  the  source  point.  It  is  difficult  to  numerically  evaluate  the  integrals  (A.4)-(A.6) 
in  this  case,  and  hence  this  singularity  is  first  removed  from  the  integrand  by  subtracting 
from  it  a  singular  term  to  which  it  is  asymptotic,  and  then  adding  back  the  analytically 
evaluated  contribution  of  the  removed  singular  term.  For  extremely  thin  wires,  better 
accuracy  of  the  numerical  integration  is  obtained  if  the  singularity  is  extracted  only  in  the 
interval  s'  €  [stt,s„],  where  su  —  max(0,  sm  —  15a)  and  sv  =  min(h,  sm  +  15a).  For 
source  points  near  observation  points,  the  elliptic  integral  behaves  as  ((17.3.26)  of  [27]) 
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(A. 15) 


The  last  term  is  extracted  from  the  kernel  for  numerical  integration,  and  its  contribution 
is  determined  by  integrating  it  analytically.  The  resulting  procedure  is  summarized  in  the 
following  section. 


A. 3  Summary  of  Computation  of  Wire  Potential  In¬ 
tegrals 


The  computation  of  potential  integrals  associated  with  a  wire  segment,  as  given  by  (A.4)- 
(A.G),  is  summarized  as  follows: 

For  an  observation  point  at  least  two  segment  lengths  away  from  the  source  segment,  the 
kernel  is  replaced  by  the  reduced  kernel  approximation,  (A. 7)  and  (A.8),  and  the  integral 


(A.16) 


is  evaluated  by  Gaussian  quadrature. 
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For  an  observation  point  less  than  two  segment  lengths  from  the  source  segment,  the 
integral  of  the  exact  kernel  is  decomposed  for  evaluation  as 

+  Hr +Qm&xw" 

+ ;  (T +  O  ™  * 


p  +  a 


I  tit  *  —  1)2,  or  3, 


where 
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-  4sm(p  —  a)  tan-1  -  a'(s' Jr  2sm)}  (. 


(A.18) 


(A. 19) 
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and  where  Rr,  sm ,  and  p  are  defined  in  (A.8),  and  (A.12),  and 

/i(a')  =  p  /,(«')  =  j  (A  -  /,(»')  =  1. 

The  quantities  su  and  are  defined  as 
su  =  max(0,  sm  —  15a), 

s„  =  min(h,  sm  +  15a). 

If  sm  is  located  outside  of  the  source  segment  then  +  J  ^  in  (A. 17)  is  replaced  by 
(0- 

All  the  integrals  in  (A.  17)  have  smooth  integrands  and  can  be  numerically  integrated 
by  Gaussian  quadrature. 
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Appendix  B 

Computation  of  Potential  Integrals 
for  Body  Basis  Functions 


The  basis  function  associated  the  nth  edge  of  a  surface  patch  model  is  given  by  (3.1), 
and  its  divergence  is  given  by  (3.2).  The  corresponding  vector  and  scalar  partial  potentials 
due  to  this  basis  function  are  given  in  terms  of  (3.14)  and  (3.15)  as  (c.f.  Fig.  3.3(b)) 
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It  is  computationally  efficient  to  simultaneously  compute  all  the  potential  integrals  asso¬ 
ciated  with  each  patch  SBq  (either  SB+  or  SB~).  The  vector  and  scalar  partial  potentials 
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Figure  B.l:  Geometrical  parameters  associated  witli  integration  of  body  and  junction  basis 
functions. 
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evaluated  at  observation  point  r  due  to  a  basis  current  directed  outward  across  the  ith 
edge  (»  =  1,2,  or  3  )  of  the  gth  triangular  patch,  SBq,  are  defined  as  (c.f.  Fig.  B.l) 
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(B.3) 


(B.4) 


If  the  current  is  directed  inward  across  the  ith  edge,  the  right  hand  sides  of  (B.3)  and  (B.4) 
are  multiplied  by  (—1). 

B.l  Treatment  of  Singularities  in  Body  Potential  In¬ 
tegrals 

The  integrands  of  Af  (r)  and  (r  )  are  singular  at  R  —  0  corresponding  to  an  observation 
point  coincident  with  a  source  point  When  r  is  not  close  to  the  source  triangle,  one 
can  easily  evaluate  (B.3)  and  (B  4i  by  numerical  quadrature  as  discussed  in  Appendix 
E.  Thus  when  r  is  close  to  or  on  the  source  triangle,  extraction  of  the  1/R  behavior  is 
necessary  before  numerical  integration  is  performed.  We  can  represent  both  of  these  cases 
by  introducing  a  closeness  parameter  a(r)  defined  as 


<r(r)  = 


0  ,  if  distance  Rc  is  larger  than  the  longest  edge  of  SBq, 
1  ,  otherwise. 


(B.5) 


where  Rc  is  the  distance  from  the  centroid  of  the  source  triangular  patch  to  the  observation 
point.  Thus  we  write  (B.3)  and  (B.4)  as 


Af(r)  =  [L  p: 
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The  integrand  of  the  first  integral  of  (B.6)  and  (B.7)  is  always  slowly  varying  and  can 
be  satisfactorily  integrated  numerically  via  the  methods  of  Appendix  E.  The  remaining 
integrals  in  both  (B.6)  and  (B.7)  can  be  integrated  analytically,  and  are  given  as  [28] 

g.  =  [  —  dS'  =  Y\  P°i  Ui  [p?  In  -  |d|  (tan-1 - T 

9‘  Jsb<  R  ^  '  *  *  R-  +  1-  11 V  (J2?)2  +  \d\R? 


(B-S) 


*  SL  TiS'  =  L^rdS' 


=  jE*  (*?)’  to  ^  |  +  W  -  l-R-  +  (p  -  p,)  9..  (B.9) 

For  numerical  evaluation  of  (B.6)  and  (B.7),  it  is  convenient  to  express  all  their  terms 
in  the  area  coordinates  of  Appendix  E.  Accordingly,  using  (E.5),  the  position  vector  p\ 
relative  to  vertex  i  is  expanded  in  terms  of  the  adjacent  edge  vectors,  which  act  as  basis 
vectors  in  the  local  coordinate  system  on  patch  SBq: 


P'i  —  Zi+lti- 1  ~ 


(B.10) 


where  the  edge  vectors  are  defined  in  Fig.  B.l.  Substituting  (B.10)  into  (B.6)  and 
using  (E.7),  we  can  express  the  vector  function  >tf  (r)  as 

A*V)  =  rz~r  [(Ab+  +  <’('•)  9*)  4-i  -  (A?-  +  o(t) g~)  <j+1] , 


(B.U) 


and 


9?  =  2A 


1  6± i 


// 
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d4±i  d£,'. 


(B.13) 


The  singular  integrals  gf  in  (B.ll)  are  related  to  the  vector  gi  of  (B.9)  as  may  be  seen  by 
substituting  (B.10)  into  (B.9)  and  comparing  to  (B.13).  We  obtain 


9i  —  9i  4-i  -  9i  4+ 1, 
where,  from  (E.10),  we  have 

^t±i 


(B.14) 


9,  =  ~ 


h. 


i±  l 


9i, 


(B.15) 


and  this  result,  together  with  (B.8)  and  (B.9),  allows  all  but  the  terms  Af*  in  (B.ll)  to 
be  evaluated. 

Similarly,  the  scalar  function  <J?f(r)  can  be  expressed  as 
^  12 


*f(  r)  = 


j4nuje  ht 


[I  +  a(r)g,  ], 


(B.16) 


with 


r i  e  —  (t(t) 

I  =  2A  /  /  - — ^  dU*  dti, 

Jo  Jo  Jl 


(B.17) 


and  g,  is  defined  as  in  (B.8). 

The  remaining  terms  in  (B.ll)  and  (B.16) — the  nonsingular  contributions  to  the  vector 
and  scalar  potentials,  Af ±  and  7,  respectively — are  computed  by  numerical  quadrature  as 
discussed  in  Appendix  E. 

B.2  Summary  of  Computation  of  Body  Potential  In¬ 
tegrals 

The  vector  potential  integral  associated  with  an  edge  as  given  by  (B.3)  is  obtained  by 
substituting  (B.8),  (B.9),  and  (B.15)  into  (B.ll)  and  is  summarized  as  follows: 
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where  <r(r)  is  defined  in  (B.5). 

The  scalar  potential  integral  associated  with  face  SBq  in  (B.4)  is  evaluated  by  substi¬ 
tuting  (B.17)  and  (B.8)  into  (B.16)  and  the  result  is  summarized  as  follows: 

„  1  9  f  /■!  p~ikR  —  rr(r\ 

>  =  -id  *  [  H  i  - ^  * +  \  ■ 

The  integrals  in  Eqs.  (B.18)  and  (B.19)  are  well-behaved  and  may  be  integrated  numerically 
via  (E.12). 


I 
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Appendix  C 

Computation  of  Potential  Integrals 
for  Junction  Basis  Functions 

The  basis  function  associated  with  current  entering  the  nth  junction  of  a  wire/body 
junction  model  is  given  by  (3.3)  and  its  divergence  is  given  by  (3.5).  The  vector  and  scalar 
potentials  due  to  a  unit  current  source  associated  with  the  nth  junction  may  be  given  in 
terms  of  (3.14)  and  (3.15)  as  (c.f.  Figs.  3.4  and  B.l) 


All  the  integrals  in  (C.l)  and  (C.2)  except  the  last  integral  in  (C.l)  are  evaluated  in 
Appendices  A  and  B.  For  the  remaining  integral,  we  first  assume  that  the  gth  junction 
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surface  patch  has  a  junction  at  vertex  t  (a  local  index)  in  the  patch.  Then  the  integral  can 
be  written  in  the  form 


C.l  Treatment  of  Singularities  in  Junction  Potential 
Integrals 

The  integrand  of  Af(r)  in  (C.3)  is  singular  when  the  source  point  is  at  junction  vertex  (i.e., 
when  p\  =  0),  and  when  an  observation  point  is  coincident  with  a  source  point  (i.e.,  when 
R  =  0).  Thus  when  r  is  close  to  the  source  triangle  or  r'  is  close  to  the  junction  vertex, 
extraction  of  the  dominant  singular  or  nearly  singular  form  of  the  integrand  is  necessary 
before  numerical  integration  is  performed.  To  discern  the  behavior  of  the  integrand  near 
these  singularities  we  expand  it  in  a  Taylor  series  about  each  singularity. 

At  p',  —  0  we  expand  terms  of  the  integrand  in  Taylor  series  in  p',  and  then  combine 
them  to  obtain  the  Taylor  series  for  the  integrand  as  a  whole.  For  example  (c.f.  Fig.  B.l), 
we  first  note  that 

R=  \Ri-  p',|  =  y/R?  +  pf  -  2ft  pV  (C.4) 

Since  for  sufficiently  small  p\,  it  is  clear  from  Fig.  B.l  that  R?  >  p'2  —  2 Hi  •  p't,  and 
hence  R  may  be  expanded  in  the  series 

R*  Ri  1  +  Y?  -  (0.5) 
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Thus  the  exponential  in  (C.3)  may  be  expanded  as 


e-jkR  ^  g-jkRi 


i  -  Jk* 


2Ri 


and  hence 


e-jkR  e-jkR, 


R  Ri 

Finally,  we  have 


1  +  (,k  +  h)  + 


Pi  e 


-jkR 


{p'i-hiY  R 


=  / 


where 


f,  = 


P,  e 


p’ 


-jkR, 


(p\ '  hi)2 


i  +  o’*  +  4-)pVHl 


(C.6) 


(C.7) 


(C.8) 


it,  7  it, 

All  the  higher  order  terms  in  the  series  vanish  as  p\  — ►  0. 

For  the  singularity  at  R  =  0,  we  expand  the  integrand  in  a  Taylor  series  in  R. 
Since  (c.f.  Fig.  B.l) 

P'i  —  Ri  —  R 


(C.9) 


(C.10) 


then 


p'i  ■  h,  =  Ri-  hi 


1  - 


and  therefore 
1 


R  hj 
Ri  hi 


,  „  R  h, 

1  +  2-^— ^  + 


(C.ll) 


Since 


{p'i -KY  {Rt-hifi  R,  h. 


e~lkR  *2it 

~7T  =  R'}k-~  + 


(C.12) 


(C.13) 
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we  have 


P'.  a  p',  1  *R  hi 

(P'i-hiP  R  (H,  •  K y  |H  ( Ri  ■  fn)R 


2 _ £i 

-fci)2 l* 


2[(p  —  Pt)  —  (p  — 

(Hi  •  fci) 


-P#)]  H  ' 

i)H 


1  f^i  .  2(P  -  Pi)  (K  ~  P'.)  ‘ 
(Hi  •  hj)2  H  (Ri-hi)R 


in  which  we  have  used  the  following  relations: 


(C.14) 


p'i  =  Ri  -  R  =  (p-  Pi)  -  (p  -  p'), 


|p  -  p'|H  in 


0  as  R  — ►  0  since  |p  —  p'\  <  R. 


All  the  higher  order  terms  in  the  series  are  bounded  and  well-behaved  as  R  — ►  0,  and  since 
R  hi  =  (H.  -  p',)  •  hi  =  (p-p<  -  p'i) -A,, 


we  define 


.  _  1  \p'i  2 (p-p,)[(p-pi)-p',]-fe. 

*  (H.-A.OUH  {Ri-hi)R 


(C.15) 


With  the  expansions  of  f  p  and  fR  in  (C.9)  and  (C.15),  then  (C.3)  can  be  rewritten  as 

A?(r)  =  [/sji  F  dS'  +  Jsj,  iS '  +  e(T)  fR  dS'  ]  ,  (C.16) 


where 


F  =  (C,7) 

and  a(r)  is  defined  in  (B.5).  The  integrand  of  the  first  integral  in  (C.16)  has  the  singularity 
at  p^  =  0  removed,  and — when  the  observation  point  is  near  the  source  triangular  patch 
SBl — has  the  1  (R  behavior  removed.  The  integrand  is  thus  bounded  and  continuous, 


< 
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and  the  integral  can  be  integrated  numerically  by  numerical  quadrature  as  described  in 
Appendix  E.  At  R  =  0,  F  is  indeterminant  and  hence,  for  numerical  purposes,  is  replaced 


by  its  limit, 


lim  F  = - ^ —  jk+  -+^'-  (2  +  jkRi)  . 

H-0  (p'.  /l,)1!  Ri  J 


(C.18) 


The  second  and  the  third  integrals  of  (C.16),  as  shown  below,  can  be  integrated  analytically. 

C.2  Efficient  Numerical  Evaluation  of  Integrals 

For  efficient  computation,  we  rewrite  the  vector  integral  equation  (C.16)  and  (C.17)  in 
terms  of  adjacent  edges  which  act  as  basis  vectors  in  a  local  coordinate  system.  Since 


p',  •  hi  =  hi  (  1  -  ), 


(0.19) 


from  (B.10)  and  (C.16-C.18),  we  have 


Ai'{r)  =  ^-Kihi 


(C.20) 


F  =  F+  /,•_ ,  -  F~  ti+u 


(C.21) 


where 


Af*  =  L  F*  JS'  +  If}  dS'  +  <7(r)  /  iS\ 

Jsj<>  JSJ1  JSJ1 


(C.22) 


1  f  ,,  P-ikR 

jt±  — - ^,i:1 _ _ _  _  _  air\  r± 


(C.23) 


»m-F±  -  + 


1  +  j  k  R{ 


(2  + jkRi)  , 


(C.24) 
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p  ■ 

/p  ~  ku  i-  e,)2  a 


—  [ 1  +  {jk  +  ~k  )^i+iti~i  ~  •  R'\ '  (c-25) 


/■ ±  =  i  f  &±i 
*  (fi<  •  fe,-)2  1  * 


2  [fc,±i  ■  (p  -  pj 

R{  •  /i,  ^*t±i 


(P  Pi)  (£<+lA-l  ^i'-l^t+l) 

i? 


(C.26) 


The  nonsingular  integrand  F*  of  (C.22)  can  be  integrated  numerically  by  Gaussian  quadra¬ 
ture  as  described  in  Appendix  E.  The  remaining  integrals  in  (C.22)  are  singular  and  may 
be  evaluated  analytically  as  follows.  From  (E.7),  we  note  the  results 

JO*  _  0/4  fl  ['-*  ii±x  Jt  Jt  _  A  //"'I  0*7\ 


L<  (i-ly  dS>  2A'l  L  (i-ly  d^t±ld^  ~  A*’ 


Similarly, 


f  _ _ dS1  =  1^.  f 

Js*  (1  -  eo2  3  ”  Js- 


£«+l£t- 1  ,0'  _  1  A  . 

/-.  f  \o  aD  ~ 


Js*  (1  -by  3  "  Js«  (l-G)2  6  ” 

where  At  is  the  area  of  the  triangular  patch  SJq .  From  (C.25-C.28),  we  now  have 

f  f±  A  1  fe-,w  (r  .  1  +  AA 

LJ-  ds  ~  7  S?  \~fiT  (6  +  ~ R~ 

(A  A-,)  (  i)  -  (A  ■/,«)  (2)  )}• 


(C.27) 


(C.28) 


(C.29) 


_J 

f  d^1  — 

2 

*«±i  •  (p  -  />,)' 

i,)2  1 

JsJo  R  Jx 

1  ’  A, 

^i±l 

From  (B.8),  (B.9),  and  (B.15),  the  singular  integrals  (C.29)  and  (C.30)  both  can  be  ob 
tained  explicitly. 


C.3  Summary  of  Computation  of  Junction  Potential 
Integrals 

The  vector  potential  integral  associated  with  a  vertex  of  a  junction  triangular  patch  as 
given  by  (C.3)  is  obtained  by  substituting  (C.22-C.23),  (C.25-C.26),  and  (C.29-C.30)  into 
(C.20)  and  is  summarized  as  follows: 


A»  =  .)• 


(C.31) 


where 


I  J/± 


2A  r  r  i-  — *ai  (*- 

'  Jo  Jo  \h>  (1-fc)2  V 


,-jkR 

~rT 


-jkR, 


Ri 

o(r) 

hi 


1  +  ijk  +  j[  )^i+lti-1  ~  6  -16+1)  •£.]) 


1 

( Ui 

2 

hi±  i 

•  (p  -  p,)' 

(Ri  ■  hty 

1  R 

J6  hi 

■ 

(p-  Pi) 

-  (6+l6-l  ~  6-1 

6+i ) 

R 

A.  f e-M  („  ,  1+jkRi 


+  cr(r) 


2\  -  l 

(Ri  •  6-1)  (  !  )-(*,•■  4+i)(  2 


f  hi±i 

2 

fr.il  •  (p  -  Pi)' 

1  ^i±  1 

5'  ii,  /l, 

h,±i 

h, 


!  i  .  a  h{+i  hi _i 

(P  -  P,)9s  +  6-1T - g,  -  6+1  T - 0, 

ft;+i  ”i— i 


(C.32) 


in  which  <r(r),  <?3,  and  <jr(  were  defined  as  (B.5),  (B.8),  and  (B.9). 
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Appendix  D 

Computation  of  Far  Fields 

We  define  a  normalized  far  electric  field  £  as 

£{r)  =  Um  r  c’kr  E\r)  (D.l) 

in  which  the  radial  phase  and  amplitude  dependence  of  the  far  field  is  removed,  leaving 
only  an  angular  dependence,  and  where,  in  the  far  field, 

E\t)  =  —ju>  [A(r)  -  A(r)  •  f  f ] .  (D.2) 

The  unit  vector  f  points  in  the  direction  of  observation  of  the  far  field,  and  is  given  by 

r  =  =  x  sin  8  cos  <f>  +  y  sin  6  sin  <f>  +  z  cos  6. 

In  the  far  field,  the  vector  potential  A(r)  reduces  to 

^(r)  =  ii  li!l  /  j(r')  eikf  r'  dS .  (D.3) 

47 rrJs 

which  is  evaluated  in  the  same  manner  as  (A.l),  (B.l)  and  (C.l)  except  that  the  kernel 
-  is  replaced  by  .  The  various  forms  that  the  integrals  take  on  are  summarized 

in  the  following  sections. 
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D.l  Far  Field  Vector  Potentials  for  Wire  Basis  Func¬ 
tions 

The  vector  potential  in  the  fair  field  due  to  a  unit  current  source  across  the  nth  node  of  a 
wire  is  given  as 

-  V~r-  4  A°(t)  e’"  T'  de‘ 

=  f  [  ~  ei"  T'  dl'  +  [  c’"  T'  dt\  (D.4) 

4n  r  |/5^+  Js?~  h™~  K 

This  integral  may  be  evaluated  in  closed  form,  but  it  is  generally  more  efficient  to  evaluate 

it  by  Gaussian  quadrature. 

D.2  Far  Field  Vector  Potentials  for  Body  Basis  Func¬ 
tions 

The  vector  potential  evaluated  in  the  far  field  point  due  to  a  basis  current  directed  outward 
across  the  ith  (t  =  1,2,  or  3  )  edge  of  a  triangular  patch  SBq  is  given  by 


Af( r)  =  f  /  £  r  dS\ 

47T  Jsbi  h, 

which,  using  (E.5)  and  (E.7),  may  be  expressed  as 

A?(-)  -  %  ~  [<4f+  4-.  -  Af-  /,+1]  , 


(D.5) 


(D.6) 


where 


Af*  =  /‘  {«,  iU,  di<. 

Jo  Jo 


(D.7) 


The  integrals  in  (D.7)  may  be  evaluated  in  closed  form,  but  it  is  generally  more  efficient 
to  evaluate  them  by  numerical  quadrature  as  described  in  Appendix  E. 
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D.3  Far  Field  Vector  Potentials  for  Junction  Basis 
Functions 

In  the  far  field,  the  vector  potential  integral  associated  with  a  vertex  of  a  junction  triangular 
patch  as  given  in  (C.3)  becomes  the  following: 


Af(r)  =  Ki  hi  —  /  — ^4—  t’kf  r'  dS'. 

*  V  Air  r  Jsj*  (p'{ .  hi )2 


(D.8) 


Note  that  there  is  still  a  singularity  at  p\  =  0,  i.e.,  when  the  source  point  is  at  the 
junction  vertex.  The  treatment  of  Appendix  C  may  be  adapted  to  remove  the  singularity 
from  the  integrand  with  the  result  that  the  final  form  of  A\'  in  the  far  field  becomes 

C  f  .J/+  t  aJ>+  a  1  rr\  r»\ 


Ai\r)  =  —  K'+  4-1  -  A?+  /,+1] , 


(D.9) 


where 


='"r7'  r  iM***  * 


0  Jo  (1  -  4)2 


=  J 1 


1  /»-*  4±1 


(1  -  6)2 


[  e>kf  r'  -  ]  d(i± ,  + 


e’kr‘Ti 


(D.10) 


The  integrals  in  the  last  line  of  (D.10)  may  be  evaluated  by  numerical  quadrature  as 
described  in  Appendix  E. 
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Appendix  E 

Local  Coordinates  and  Numerical 
Integration  on  Triangular  Patches 


In  this  Appendix,  we  introduce  a  convenient  set  of  local  coordinates  defined  on  triangular 
patches.  This  set  of  coordinates  is  useful  for  representing  both  scalar  and  vector  functions 
on  triangles  and  for  performing  numerical  integration  over  the  patches. 


E.l  Normalized  Area  Coordinate  System 


We  introduce  the  so-called  normalized  area  coordinates  £,  (i  =  1,2,3),  defined  as 


ii  = 


A, 

A  ’ 


(E.l) 


■where  the  area  A,  is  the  area  of  the  triangle  formed  by  point  r  and  edge  i  (Fig.  E.l),  and 
A  is  the  total  area  of  triangle  Sq.  Hence 

£(.  =  1.  (E.2) 

1=1 

Cleaxly,  only  two  of  the  coordinates  axe  independent,  and  the  third  can  always  be  elimi¬ 
nated  via  (E.2).  The  transformation  between  the  global  and  local  coordinate  systems  is 
summarized  by  the  relation 


r  =  &+lri+J  +  ir«— i  +  («»\ 


(E.3) 


84 


1 


where  the  position  vector  r  and  the  vertex  vectors  rj,  1*2,  and  f3  are  given  in  the  global 
coordinate  system.  Using  (E.2)  to  express  the  position  vector  in  terms  of  the  two  indepen¬ 
dent  coordinates  £,±j  and  and  noting  that  r,±i  —  r,Ti  =  and  r^i  —  r,-  =  we 

have 


*•  =  (!  —  1  —  &)r*Tl  +  &±lri±l  + 


=  _  &*.±l)  +  *»*!• 


(E.4) 


A  frequently  used  result  is  that 


P'i 


r'  -  r, 

(i+ir,+i  +  +  (fr  -  l)r, 


=  £«+ll\+l  +  ~  (&+1  +  £«+l)ri 

=  £*+l^i-l  —  {i-l^i+1-  (E.5) 

From  (E.4)  we  also  have 

dS,  =  x  dli 

—  |^«  ^  t-i-kl  |  1  d£i 

=  2A  dUi  dii.  (E.6) 

From  (E.5)  and  (E.6),  it  follows  that  surface  integrals  over  an  arbitrary  scalar  function 
f(p\)  on  Sq  can  be  expressed  in  terms  of  local  coordinates  as 

[  f(p't)dS'  =  2 A  [l  r1,  ,  -  (i-iii+i)  dUi  dti-  (E.7) 

JS«  JO  JO 


E.2  Expansion  of  Tangent  Vectors  on  a  Triangle  in 
Terms  of  Edge  Vectors 


Any  tangent  vector  in  S?  can  be  expressed  in  terms  of  a  pair  of  triangle  edge  vectors,  which 
act  as  basis  vectors  in  the  local  coordinate  system.  For  an  arbitrary  tangent  vector  /,  for 
example,  we  wish  to  write  it  as 
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i.e.,  we  wish  to  determine  its  expansion  coefficients  f+  and  /  in  terms  of  (non-orthogonal) 
basis  vectors  l,-i  and  l,+i .  This  is  easily  done  once  we  notice  that  (c.f.  Fig.  E.2) 


hi  ti  =  0,  (E.9) 

(i.e.  l{  and  hi  are  biorthogonat)  where  hi  —  hihi.  Hence,  from  (E.8)  and  (E.9),  and  since 
hi± i  ••^1^1  =  T^{±i  » 

we  have 


Thus 


hi±i  f  =  ±  /*  hi± i  •  =  -  Z1  /i-±1. 


/*  =  -  ^ 


^«±i 

and  finally,  from  (E.8), 


(E.10) 


f  ~  f  b-lhj+l  +  A+lfe,  -  y. 


(E.ll) 


\  K  +  \  hi_  i 

where  the  quantity  in  parentheses  is  the  identity  dyad  in  local  coordinates,  multiplication 
by  which  expresses  a  tangent  vector  in  terms  of  the  triangle’s  edge  vectors.  Since  the 
vectors  li±x  in  (E.8)  are  constant,  a  frequent  use  of  this  result  is  to  reduce  an  integral  over 
an  arbitrary  vector  f  to  two  scalar  integrals  of  the  form  of  (E.7)  over  the  components  f±. 


E.3  Numerical  Integration  over  Triangular  Patches 

Numerical  integration  over  a  triangle  can  be  accomplished  in  terms  of  the  area  coordinates 


(E.12) 


/'  n  /te. fi«)  <<{•  <%*i  *  £ 

«/o  */o  J~1 

where  iV  is  the  number  of  evaluation  points  of  the  quadrature.  Appropriate  weight  coeffi¬ 
cients  u>j  and  coordinates  ((&);,(&±i)j)  of  points  at  which  the  function  must  be  evaluated 
for  Gaussian  schemes  are  given  in  [29].  Only  one-point,  three-point,  and  seven-point 
schemes  (depending  on  the  ratio  of  distance  between  observation  point  and  centroid  to  the 
maximum  edge  length  of  the  source  triangular  patch)  are  used  in  the  code  JUNCTION, 
and  the  corresponding  weight  coefficients  and  corresponding  coordinates  axe  presented  in 
Tables  E.1-E.3. 
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u 

[• 


Y' 


Table  E.l:  Weight  Coefficients  and  Local  Coordinates  for  One-Point  Gaussian  Quadrature 
on  a  Triangle. 


j  tVj 

(fc)i 

(fc±i)i 

0.333333333333333 

0.333333333333333 

Table  E.2:  Weight  Coefficients  and  Local  Coordinates  for  Three-Point  Gaussian  Quadra¬ 
ture  on  a  Triangle. 


j 

Wj 

(fc)i 

1 

0.166666666666667 

0.666666666666667 

0.166666666666667 

2 

0.166666666666667 

0.166666666666667 

0.666666666666667 

3 

0.166666666666667 

0.166666666666667 

0.166666666666667 

Table  E.3:  Weight  Coefficients  and  Local  Coordinates  for  Seven-Point  Gaussian  Quadra¬ 
ture  on  a  Triangle. 


j 

w} 

(fc); 

1 

0.1125 

0.333333333333333 

0.333333333333333 

2 

0.062969590272413 

0.797426985353087 

0.101286507323456 

3 

0.062969590272413 

0.101286507323456 

0.797426985353087 

4 

0.062969590272413 

0.101286507323456 

0.101286507323456 

5 

0.066197076394253 

0.470142064105115 

0.470142064105115 

6 

0.066197076394253 

0.470142064105115 

0.059715871789770 

7 

0.066197076394253 

0.059715871789770 

0.470142064105115 
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